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Abstract 

We describe a method to approximate the matrix elements of a real 
power a of a positive (for a £ R) or non- negative (for a > 0), infinite, 
bounded, sparse and Hermitian matrix W . The approximation uses only 
a finite part of the matrix W . 
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1 Introduction 

The motivation for this paper has been the research, carried out by Jerzy Ki- 
jowski [6l [5J |4j and, for the last year, the author, on the quantum field theory 
on a lattice. This physical theory arises from the discretization of the classical 
field Lagrangian and its subsequent quantization. As the number of degrees 
of freedom is infinite, Stone-von Neumann's theorem [8, 9|, which proves the 
uniqueness (up to a unitary transformation) of the representation of the algebra 
of field operators, does not apply. The problem of finding such a representa- 
tion has been customarily solved by resorting to global symmetries, such as the 
translational invariance. Our research has been concerned, however, with the 
quantum field theory on the background of non-homogeneous gravitational field. 
This eliminates global symmetries and renders traditional ways of finding the 
physical representation obsolete. 

An attempt was made to find a representation of the field operators on a 
Hilbert space via a limit of the representations on finite lattices, which are, 
as proven by von Neumann, uniquely defined up to a unitary transformation. 
When applied to a quadratic Hamiltonian on a c?-dimensional rectangular lattice: 

\ Y + Y Y W ™™(i> m <j> n , 
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(where n is a vector of integer indices numbering the lattice sites, 7r n is the 
(real) momentum operator for the site n, <j) n is the field operator, also for the 
site ra, and W an infinite, positive, real and symmetric matrix), it turns out that 
the existence of the limit of representations is connected with the existence of 
the following limit: for each finite lattice T C Z d , we define a finite matrix Wr, 
which is a finite part of W corresponding to the lattice T. In order to make the 
problem well-posed physically, we need to impose some boundary conditions on 
the elements of Wr corresponding to the edges of T. To ensure the existence of 
the groundstate, we choose non-negative boundary conditions. It turns out that 
the existence of the limit of finite-dimensional representations is guaranteed by 
the convergence of the sequence of matrix elements (W T ) m n with increasing 
lattices T. Thus, the problem is a subset of a more general problem: for a 
real power of an infinite, positive and bounded Hermitian matrix W, we try to 
approximate its matrix elements with the corresponding elements of its finite 
parts, with possible boundary conditions added. This is the topic of this paper. 

2 Main result 

An infinite matrix W is also a linear operator on the Hilbert space £ 2 of square- 
summable infinite sequences. The matrix elements W mn are equal to the scalar 
product {We n \e m ), where the e„ are the canonical basis vectors of I 2 . Con- 
versely, any bounded linear operator on I 2 corresponds to a certain infinite 
matrix. 

The infinite matrix W a is, like W , a positive, bounded, linear and Hermitian 
operator on I 2 which, when acting on any eigenvector v of W with an eigenvalue 
A, fulfills the equation W a v — X a v. For a general theory of fractional powers of 
operators, see [2J. 

The core result of this paper is 

Theorem 1 Consider an infinite Hermitian matrix W with the following prop- 
erties: 

1. \\W\\ <oo, 

2. each row of W has at most k non-zero elements, 

3. inf a(W) > 0, where a(W) is the spectrum ofW as an operator on I 2 . 

For a 6 K (for inf a{W) > 0) or a > (for inf a(W) = 0), we have for a 
fixed pair of indices m,n £ Z 

lim (W$ Q ) mn = W% n (pointwise). (1) 

P,Q— >oo ^ 

where Wpq is a (P + Q + 1) x (P + Q + 1) matrix with indices in the range 
[— P, Q] which is defined as follows: 

(W PQ ) mn = W mn , -P <m<QV -P <n<Q (2) 
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{Wp Q ) mn = W mn + D mn , m,ne {-P,Q} . (3) 

The Hermitian matrix D fulfills the following conditions: for any P,Q>0, we 
have 

W PQ > inftr(W) (4) 
(i.e. the lowest eigenvalue ofWpQ is greater or equal to inf a{W)) and 

\\W PQ \\ < \\W\\+ d, (5) 

for d>0. 

3 Physical example 

Consider an infinite matrix M which generates the following quadratic form for 

x e I 2 : 

oo oo 

x^Mx = a ^2 \x n \ 2 +b ^2 \ x n - x n-i\ 2 , a>0Ai)>0. 

n=— oo n— — oo 

It arises from the homogeneous theory of scalar field on a one-dimensional lat- 
tice. The finite matrix Mpq is chosen in such a way that it generates the 
following quadratic form for x G C p+< 3 +1 : 

Q Q 
x^Mpgx^a ^ \ x n\ 2 +b ^2 l x « ~ x n-i\ 2 + b \x Q - x_ P \ 2 , 

x=-P n=l-P 

an example of the so-called "periodic boundary conditions" in physics. It turns 
out that for — P < m,n < Q and aelwc have 

{MpQ)mn - ^ P + Q + l C ° S [p+Q+l) ■ 
On the other hand, 

M rnn = [ ( a + b s ' m2 ^k) a cos[27rfc(™ - n)]dk . 
Jo 

It is clear that 

Hm(M£ Q ) mn = M«„ . 
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4 Proof of Theorem [T] 

For a = 0, the theorem is trivial. 

Since every row of W has at most k non-zero elements, any matrix element 
of any positive and integer power of W can be calculated as a finite sum of finite 
products of the matrix elements of W: 



mi Gr(m) 

E 



E E W mmi W mim2 W m2n (6) 



where 



r(m) := {n|W m „ ^ 0} . 

Because of ([2]) and J6|), for any indices m,n and for every non- negative and 
integer power j one can choose such finite Pj and Qj that 



yp> P] .Q>Q 3 Wl n = {Wpq) 



ftlli 



This proves for cieN. 

Because of ©, ©, ((U) and ([5]), we have the following bounds for the spec- 
trum of Wpq as an operator on C pq , the space of all complex finite sequences 
indexed from —P to Q: 

o-pq{W pq ) C [Ma{W),\\W\\+d] . 

We make the distinction between o~pq (spectrum of an operator on C-pq) and 
a (spectrum of an operator on £ 2 ), because Wpq, as an operator on £ 2 , has 
infinitely many eigenvectors with the zero eigenvalue: 

W PQ e n = 0, n < -P V n > Q , 

which means that 

a(W P Q) = ap Q (Wp Q )U{0} . 
The lower and upper bound will be denoted by c and w, respectively: 

c := inf a{W) , c > 
w := \\W\\ + d , c < w < oo . 

Since \\W\\ < w, W a can be expressed [7] as the following power series: 

where / is the infinite identity matrix, I mn — S m ^ n . We know that \\W — wl\\ < 
w — c < w. Hence, \\W — wl\\ /w < 1 and the series ([7]) is weakly convergent (in 
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the case of \\W — wl\\ /w = 1 the series is weakly convergent for a > only). 
Weak convergence of (J7J implies the convergence of the following series: 



(W a ) r 



a\ (W - wiy m 
i, 



(8) 



where {W a ) m n is defined as a matrix element of the linear operator W a between 
two canonical basis vector of i 2 , e m and e n . 

Since every row of W has at most k non-zero elements, any matrix element 
of any positive and integer power of W — wl can be calculated as a finite sum 
of finite products of the matrix elements of W — wl: 



2 

mn 



3 

mn 



(W - wI) mmi (W - wl) 

m\ Gr'(m) 

Y, E (W ~ wl) mmi (W - wl) mim2 (W - wl), 

m±Gr' (m) m2£r'(mi) 



(9) 



(W - wl) 
(W - wl) 
(W - wlf n 

where 

r'(m) := {n\(W - wl) mn ^ 0} . 
Because of © and © , we can define a non- negative integer jpq such that 

i < jpq (w - wiy mn = (w PQ - wi PQ y mn , (10) 

where Ipq is the (P + Q + 1) x (P + Q + l) identity matrix, and 



lim 7 on = oo . 

P,Q-oo w 



(11) 



(Wpg) TOn , for Wpq being an operator on Cpq, can be expressed as a series 
similar to (El): 



m Q ) 



W \ j J wi 



)l 

'inn 



We split this scries with the help of ([TO 

jpq— i 



J=0 



inn 



W 



E 



cA (Wpq - wi PQ y r 
j) 



wJ 



Using (0) and (fl2|l . we get 



(12) 



3=3PC 



3=3PQ 



3 

inn 



jj 
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The difference can be estimated as follows: 
\(W£ Q ) mn -(W a ) mn \ = 



y, /g\ (W P Q - wI PQ y mn y, / g\ (W - wl) 



< w° 



3=3 P 
oo 



E 

3=3 PQ 



\{W P q - wI PQ ) 



1 

mil 



3=3PQ 



E 

i=jpQ 



3 

run 



(W - wl)i 



(13) 

It is a well-known fact that, for every Hermitian operator A on a Hilbert space 
TL which has a countable orthonormal basis {e„}, every matrix element of A 
fulfills the inequality 

\A mn \ < \\A\\ . 
Indeed, using the Cauchy-Schwarz inequality we obtain 

\A mn \ = |(Ae n |e m )| < ||e m || \\Ae n \\ < \\A\\ \\e m \\ \\e n \\ = \\A\\ . 
Since £ 2 has a countable orthonormal basis and W — vol is Hermitian, we have 

\(w-wiy mn \ <( W - C y . 

Ditto 

\(w PQ - wi PQ y mn \ <(w- c y , 

since Wpq — wIpq is a Hermitian matrix. We insert these bounds into (fl~3|) and 
obtain 



\{W? Q )mn ~ (W a ) mn \ < 2w a J2 Q 



3=3PQ 



The series 



E 

3=0 



(14) 



(15) 



(i.e. the same as above with jpQ set to 0) converges (as with (0, for c = it 
converges only for a > 0). For negative a and c > 0, 



E 

3=0 



E 

3=0 



a\ I c — w 
3 



w 



{c/w) a , a < 



(16) 



For positive and non-integer a, we have ([a] denotes the integer part of a): 



E 

3=0 
W+l 



£ ^ (ro - e)i + ( - 1)M(e -,) j +( _ i)M+1(e/rori q>oa ^ 



Z 



(17) 



G 



Convergence of (fTSj) , proven in ([To]) and p7|) , shows that the sum of the elements 
of the series (| 13[) from j = jpQ to oo goes to zero for jpQ — * oo: 



lim 2w a V [" 

3=3PQ 

Using dn}, © and JTSJ, we calculate 



w — c 



(18) 



lim |(W%) m „-(ir%J < lim 2 W Q V 



J=Jpq 



lim 2tO Q V 
3=3PQ 



which proves |T]). 



w — c 



w — c 



if 



, 



5 Accuracy and applicability of the approxima- 
tion 

Inequality (fTi)) and results ([l7j|) and (|17p can be used to estimate the accuracy of 
our approximation. For a given unknown matrix element W^ n , we are interested 
in how fast the right side of (|14|1 diminishes with increasing P and Q. This 
involves two problems: 

1 . the relation between P, Q and jpQ , which is determined by the structure 
of W near the chosen matrix element, and 

2. the speed of convergence of the series (JT5J) , which appears in (fT4|) . 

Since the series (|15[) has only positive elements, its convergence is determined 
by how fast they vanish with increasing j. This in turn is determined by the 
power factor (1 —c/w)^ (for c > 0) and the binomial coefficient (") . The smaller 
is 1 — c/vu, the faster will (fT5| converge. The dependence of (") on j is such 
that (I15|) will converge more slowly when a < — 1, than in the opposite case. 

The relation between P, Q and jpQ (i.e. how fast jpQ diverges when P, Q — > 
oo) cannot be determined without additional assumptions about W . It is ob- 
vious that jpQ = 1 (the minimal value it can take) for — P > min(m, n) or 
Q < max(m, n), where m,n are indices of the matrix element that we want to 
compute. 

Suppose, for example, that W is fc-diagonal for k = 21 + 1, so that W mn = 
for |m — n| > I and is non-zero otherwise. We then have jpQ = l+[min(min(m, n)- 
P — 1, Q — max(m, n) — for — P < min(m, n) — 1 and Q > max(m, n) + 1. 
Asymptotically, j'pq grows linearly with min(P, Q). 

Since W is sparse, positive integer powers may be easily computed without 
resorting to any approximations, using (|6|). 



7 



Our approximation can be used to approximate locally the solutions to infi- 
nite matrix equations Wx = f (i.e. approximate the solutions of linear operator 
equations, see [DI3j), for sequences / which are non-zero only for a finite subset 
of indices. Additionally, the approximation might be applied also to calculate 
numerically the matrix element of some power of a finite, but very large matrix 
which fulfills the conditions of our theorem. 

6 Summary and generalization 

We have shown that one can use a finite part of the infinite, non-negative, Her- 
mitian, sparse bounded matrix W to approximate a chosen matrix element of a 
power W a for any real or positive a (for inf W > or inf W = 0, respectively). 
The accuracy of the approximation depends on the bounds of the matrix spec- 
trum, the value of a and on the structure of W near the chosen element. Possible 
applications of the approximation have been discussed briefly. 

It is apparent that the result presented in this paper can be easily generalized 
to the case of not only power functions, but also all functions which arc finite and 
possess a power series expansion absolutely convergent on the whole spectrum 
of W. 
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